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^ . We report the results of spectroscopic mapping observations carried out to- 

ward the Herbig-Haro objects HH7-11 and HH54 over the 5.2 — 37/im region 
using the Infrared Spectrograph of the Spitzer Space Telescope. These observa- 
tions have led to the detection and mapping of the S(0) - S(7) pure rotational 
CN ! lines of molecular hydrogen, together with emissions in fine structure transitions 

^ ! of Ne^, Si^, S, and Fe"*". The H2 rotational emissions indicate the presence of 

\o 

O . 

. tent with the expected temperature behind nondissociative shocks of velocity 

~ 10 — 20kms~^ - while the fine structure emissions originate in faster shocks 
of velocity ~ 35 — 90kms~^ that are dissociative and ionizing. Maps of the 
H2 line ratios reveal little spatial variation in the typical admixture of gas tem- 
peratures in the mapped regions, but show that the H2 ortho-to-para ratio is 
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j§ ! quite variable, typically falling substantially below the equilibrium value of 3 

attained at the measured gas temperatures. The non-equihbrium ortho-to-para 
^ I ratios are characteristic of temperatures as low as ~ 50 K, and are a remnant 

^1 of an earlier epoch, before the gas temperature was elevated by the passage of 

a shock. Correlations between the gas temperature and H2 ortho-to-para ratio 
show that ortho-to-para ratios < 0.8 are attained only at gas temperatures be- 
low ~ 900 K; this behavior is consistent with theoretical models in which the 
conversion of para- to ortho-H2 behind the shock is driven by reactive collisions 
with atomic hydrogen, a process which possesses a substantial activation energy 
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barrier {EA/k ~ 4000 K) and is therefore very inefficient at low temperature. 
The lowest observed ortho-to-para ratios of only ~ 0.25 suggest that the shocks 
in HH54 and HH7 are propagating into cold clouds of temperature ^ 50 K in 
which the H2 ortho-to-para ratio is close to equilibrium. 

Subject headings: ISM: Molecules — ISM: Abundances — ISM: Clouds - molec- 
ular processes - shock waves 

1. Introduction 

Stars have a profound effect upon the interstellar medium (ISM) that surrounds them. 
Stellar radiation heats - and sometimes photoionizes and photodissociates - the interstellar 
gas, while supersonic outflows associated with stellar birth (in protostellar outflows) and 
death (in supernova remnants) send shock waves propagating through the ISM. In regions 
of high mass star-formation - such as the Orion Molecular Cloud, for example - many of 
these processes operate simultaneously, making their separate effects hard to disentangle and 
elucidate. In Herbig-Haro objects and supernova remnants, by contrast, the effects of shock 
waves can be studied in relative isolation in an environment where the stellar radiation field 
is not greatly enhanced. 

Over more than three decades, the physics and chemistry of interstellar shock waves 
have been the subject of intensive study, both theoretical and observational, the results 
of which have been summarized in several review articles (e.g. Walmsley & Pineau Des 
Forets 2005; Draine & McKee 1993). When shock waves propagate within molecular clouds, 
the shock velocity and magnetic field determine the fate of the gas. Fast shocks (with 
propagation velocities Vg ^ 40kms~^ given interstellar magnetic fields of typical strength) 
are dissociative: they destroy molecules and may ionize the resultant atoms. Slow shocks 
{vg ^ 40kms'"^), by contrast, are non-dissociative: they heat and compress the molecular gas 
through which they pass, and may modify its chemical composition, but molecules survive 
their passage. The chemical changes wrought by non-dissociative shocks typically result 
from elevated gas temperatures or the phenomenon of ion-neutral drift (crucially important 
in non-dissociative shocks), which can drive chemical processes that are negligibly slow in 
cold clouds because they are endothermic or possess activation energy barriers. The result of 
these processes depends critically upon the timescale for key chemical reactions relative to the 
period of time for which shock elevates the temperature and ion-neutral drift velocities. In a 
previous study of HH54 carried out with the Infrared Space Observatory (ISO), Neufeld et al. 
(1998; hereafter NMH98) argued that the H2 ortho-to-para ratio, observed to he significantly 
below its equilibrium value (~ 3) in the warm {T ~ 650 K) shocked gas in HH54, provided 
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a valuable probe of the relevant timescales. 

The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope is a powerful 
tool for the study of shocked interstellar gas. With a wavelength coverage of 5.2 - 37 /xm, 
Spitzer/IRS provides access to the S(0) through S(7) pure rotational lines of H2, as well as 
fine structure emissions of several atoms and ions produced in fast, dissociative shocks: Ne"*", 
Si"*", S, and Fe"*". The great sensitivity of Spitzer/IRS, together with the multiplex advantage 
inherent in long slit spectroscopy, allows large regions to be mapped rapidly. 

In this paper we report the results of spectral line mapping in HH54 and HH7 — 11, two 
well-studied Herbig-Haro objects known to contain warm molecular gas. HH 7-11 is a chain 
of Herbig-Haro objects in the Perseus molecular cloud, associated with NGC 1333 (Snell & 
Edwards 1981). These objects lie in the blue-shifted part of a CO outflow that originates 
near the protostar SVS 13. HH 7-11 has been studied extensively through observations 
of optical emission lines (Hartigan, Curiel, & Raymond 1989), near-IR emission lines of 
[Fe II] (Stapelfeldt et al 1991; Gredel 1996), H2 vibrational emissions (Noriega-Crespo et 
al. 2002, Khanzadyan et al. 2003), far- infrared atomic and molecular lines (Molinari et al. 
2000); and through high resolution studies of CO emission (Bachiller et al. 2000). One 
key result to emerge from previous studies is that the hot H2 emission exhibits a classic 
bow shaped morphology in HH 7, with an emission peak offset from the peak emission 
of atomic flne structure lines. Observations of the HH7-11 flow have proven valuable in 
constraining detailed models for shock excited gas (e.g. Hartigan, Raymond, & Hartmann 
1987; Smith, Khanzadyan, & Davis 2003). The distance to HH7-11 is a matter of debate; 
we adopt the value of 250 pc used by Enoch et al. (2006; see discussion in that paper which 
refers to Cernis 1993; Cernis & Straizys 2003; Belikov et al. 2002). HH54 is located in the 
Chamaeleon II molecular cloud at an assumed distance of 200 pc (Hughes & Hartigan 1992), 
and is known to have a large column of warm (> 2000 K) shocked molecular and atomic gas 
(Schwarz & Dopita 1980; Gredel 1994). While it is less well studied than HH7-11, HH54 is 
notable because of an unusually large abundance of water vapor (Liseau et al. 1996) and 
H2 ortho-to-para ratios that are below the equilibrium value expected at the observed gas 
temperature. 

Our observations and data reduction methods are discussed in §2 and Appendix A, and 
the results are presented in §3. The implications of these results for our understanding of 
interstellar shock waves are discussed in §4. A brief summary follows in §5. 
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2. Observations and data reduction 

One field in HH54, and two fields in HH7-11 - one centered on HH7 and one on the 
protostar SVS13 - were observed using the IRS as part of the IRAC Guaranteed Time 
program. The Short-Low (SL), Short-High (SH) and Long-High (LH) modules were used 
to provide complete spectral coverage over the wavelength range available to IRS. Maps 
covering fields of size ~ 1' x 1' were obtained by stepping the slit perpendicular, and - in 
the case of SH and LH - parallel to its long axis. The observational details are provided in 
Table 1. 

The data were processed at the Spitzer Science Center (SSC), using version 12 of the 
processing pipehne, and then reduced further using the SMART software package (Higdon et 
al. 2004), supplemented by additional routines that we have developed. The key capabilities 
provided by these additional routines were the removal of bad pixels in the LH and SH data, 
the calibration of fiuxes obtained for extended sources, the extraction of individual spectra 
at each position sampled by the IRS slit, and the creation of spectral line maps from the set 
of extracted spectra. Details of our data reduction method, as well as the various tests we 
performed to validate the additional routines that we developed, are discussed in Appendix 
A. 



3. Results 

3.1. Spectral line maps 

Widespread rotational emissions from molecular hydrogen, along with fine structure 
emissions from atomic sulfur and several singly-ionized species, were detected in each of the 
two sources we mapped. In Figures 1 and 2 we present maps of the H2 line emissions detected 
in HH54. The straight lines demark the regions within which each transition was mapped. 
The S(l) through S(7) lines all show morphological similarities, with most maps showing 
local maxima close to chimps E and K and to clump C (in the designation of Sandell et 
al. 1987). However, real differences in the spatial distribution of the various transitions are 
apparent and provide key information about the spatial variation of the H2 excitation and 
ortho-to-para ratio; these variations are discussed in detail in §4.2 and §4.4 below. 

Emissions in several fine structure transitions have also been detected toward HH54. 
Figure 3 shows the distribution of the four fine structure emission lines that are strong 
enough to map: [Nell] 12.8/xm, [Fell] 26/xm, [SI] 25/xm, and [Sill] 35/xm. Offsets between 
the fine structure emissions and the H2 emissions are clearly evident in Figure 7 (left panel) , 
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which shows contour plots that compare the distributions of the [Nell] and H2 S(3) emissions, 
suggesting a different origin for the H2 and fine structure emissions, discussed further in §4.1 

Figures 4 and 5 show maps of the H2 hne emissions detected in HH7-11. As before, 
the straight hnes demark the regions within which each transition was mapped; in this case, 
two overlapping rectangular fields were mapped in the SL and SH modules and one field in 
the LH module. The dashed circles centered at offset (0,0) encircle the regions close to the 
bright contimmm source SVS13 where the line intensity estimates are unreliable. The S(l), 
S(2), S(3) and S(5) line maps are remarkably extended, with diffuse emission covering (and 
extending up to 20" on either side of) the line joining HH7 to SVS13. A detailed discussion 
of the differences in the spatial distribution of the various transitions - and of the inferred 
variations in the H2 excitation and ortho-to-para ratio - is again deferred to §4.2 and §4.4. 

As in HH54, emissions in several fine structure transitions have also been detected in 
HH7-11. Figure 6 shows maps of the [Fell] 26/im, [SI] 25/xm, and [Sill] 35/im line emissions 
detected with the LH module toward HH7. Unlike in HH54, the [Nell] 12.8 /xm line is not 
detected. Figure 7 (right panel) shows that the fine structure emissions, as traced by [Sill] in 
this case, are again offset from the H2 emissions, peaking several arcseconds closer to SVS13. 

3.2. Average spectra 

We have averaged the observed spectra over several regions within the HH54 and HH7-11 
fields, thereby obtaining substantial improvements in the signal-to-noise ratio. The resultant 
spectra, plotted in Figure 8, are the averages of all spectra within a series of circular regions, 
the contributing spectra being weighted by a Gaussian taper (HPBW of 15") from the center 
of each synthesized beam. The position of each synthesized beam is demarked by the dotted 
circles in Figure 7, and the central positions are given in the caption to Figure 8. 

In HH54, we obtained average spectra for circular regions centered on the peak of the 
[Nell] emission and two peaks of the H2 emission. We denote the three circular beams 
HH54FS (the peak of the [Nell] emission), HH54C (the northern of the two H2 peaks, being 
roughly coincident with clump C identified by Sandell et al. 1987), and HH54E+K (the 
southern H2 peak, roughly coincident with clumps E and K). In IIII7-11, we obtained the 
average spectrum for a circular region centered on the HH7 bow shock. 

The spectra plotted in Figure 8 all show very weak continua, leading to large line-to- 
continuum ratios even in the low- resolution module (Short-Low, left panel, with A/ A A of 
only ~60). A fairly rich spectrum of H2 rotational lines and low-excitation fine structure 
lines is apparent in each case. 
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A total of eight H2 rotational lines and eight fine structure fines of Ne+, Si+, S, and Fe+ 
(5 transitions) have been detected at one or more such region. Table 2 lists the detected 
lines together with the mean frequency-integrated intensities for each circular region. As 
discussed in §4.1 below, the H2 line ratios and intensities are typical of slow nondissociativc 
shocks, while the fine structure line intensities are in good agreement with those expected 
from fast dissociative shocks. 

No fewer than five [Fell] fine structure lines are clearly detected toward position HH54FS. 
The observed spectra are shown in Figure 9, and the transitions detected are marked on the 
Grotrian diagram in Figure 10. The strength of the higher excitation transitions originating 
in the term, relative to the lower transitions within the term, is expected to be an 
increasing function of both temperature and density. As discussed in §4.1, the [Fell] line 
ratios provide a valuable probe of the physical conditions in the emitting region. 

Upper limits are also obtained for two rotational transitions of water vapor that lie 
within the LH band and are discussed in §4.3 below. In addition, the R(3) and R(4) rotational 
lines of HD have been tentatively detected and discussed elsewhere (Neufeld et al. 2006). 

4. Discussion 

4.1. Fast and slow shocks 

In both sources that we observed, clear offsets between the fine-structure emissions 
and the molecular hydrogen emissions are evident. Similar offsets have also been observed 
toward Herbig-Haro objects in the GGD37 (Cepheus A West) region (Raines, 2000). The 
fine structure emissions in HH7 are offset by ~ 5" from the bow shock, corresponding to a 
projected distance of ~ 2 x 10^^ cm for an assumed distance to the source of 250 pc. Such 
an offset has been observed previously in HH7 (Staplefeldt et al. 1991), where near-infrared 
emissions from [Fell] and H2 have been mapped at subarcsecond resolution by Noriega- 
Crespo et al. (2002) using the Hubble Space Telescope. Noriega- Crespo et al. interpreted the 
[Fell] emissions as emerging from a "Mach disk" , in which the colfimated supersonic outfiow 
is decelerated and heated in a fast shock. In this picture, the fast shock is offset from a 
slower non-dissociative bow shock that is driven into the ambient molecular medium, the 
latter giving rise to the molecular hydrogen emissions observed from HH7. Our mid-infrared 
observations indicate that the [Sill] and [SI] fine structure emissions - like the [Fell] emissions 
- also emerge from a region offset from the bow shock. The morphology of the [Sill] map is 
most similar to the [Fell] map. These observations support the two-shock picture in which 
the [Fell], [SI] and [Sill] fine structure emissions emerge from fast shock. 
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Our observations of HH54 indicate a similar offset (~ 10", corresponding to a projected 
distance of ~ 3 x 10^^ cm for a source distance of 200 pc) between the fine structure and 
molecular hydrogen emissions. As in HH7, the ccntroid of the fine structure emissions lies 
closer to the presumed source of the outflow than the H2 emissions, consistent with the 
two-shock picture in which the fine structure emissions trace a fast shock that decelerates 
the collimated outflow. 

We have compared the observed fine structure line intensities listed in Table 2 with pre- 
dictions of the fast shock model of Hollenbach & McKee (1989, hereafter HM89). Because 
the beam filling factor of shocked material is not known, line ratios provide a more valuable 
diagnostic than absolute intensities. Although the predictions of the HM89 shock models 
are undoubtedly somewhat uncertain - given substantial uncertainties in the coUisional ex- 
citation rates for the various transitions and in the gas-phase abundances of the various 
elements - the predicted [Nell] 12.8//m/[FeII] 26.0 //m fine ratio is a strongly increasing 
function of the shock velocity, values ^ 60kms~^ being needed to ionize Ne significantly. 
The [Nell] 12.8 /im/ [Fell] 26.0 /xm line ratio of ~ 0.6 observed in HH54 requires a shock 
velocity ~ 80kms^^. In HH7, by contrast, the absence of detectable [Nell] 12.8 emission 
places an upper limit Vg < 70kms^^, whilst our measurement of the [Fell] 17.9 /xm transition 
requires a shock velocity^ of at least 35kms~^. 

The HM89 shock models predict that the [Fell] and [Sill] line intensities show a very 
similar dependence upon the shock velocity, whilst the [SI] /[Fell] and [SI] /[Sill] line ratios 
are a decreasing function of shock velocity. This prediction provides a natural explanation 
for the fact that the [Sill] map is most similar to the [Fell] map: insofar as the fine structure 
emissions may arise from a superposition of shocks at varying velocities, the [SI] emissions 
would be expected to show a spatial distribution that is somewhat different from [Sill] and 
[Fell]. 

An independent check of the two-shock picture is provided by our detection of multiple 
[Fell] lines. The mimimum temperature derived from the [Fell] line ratios (~ 5000 K) is 
larger than the maximum temperature at which H2 can survive for any appreciable time 
period - and is a factor ~ 5 larger than that inferred from the H2 line ratios (see §4.2 below) 
- lending strong support to the theory that both dissociative and nondissociative shocks are 
present in the source. 



"'^Considerably smaller shock velocities would be obtained in a comparison with the model predictions of 
Hartigan et al. (1989), which apply to shocks in low density clouds that are assumed to be ionized prior to 
the passage of the shock. 
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4.2. The temperature and column density of shocked molecular hydrogen 

The H2 S(0) - S(7) emissions provide a valuable probe of the temperature distribution 
in the shocked gas in which H2 is the dominant constituent. The observed line intensities for 
these H2 quadrupole transitions provide a direct measurement of the column density, iV( J^) 
in each upper state, J^, because all the hues involved are optically-thin. As in NMH98, we 
applied extinction corrections based upon the interstellar extinction curves of Weingartner 
& Draine (2001; "Milky Way, Ry = 3.1"). For HH54, we adopted Gredel's (1994) estimate 
of 0.3 for E{J — K), obtained from observations of the [Fell] 1.257 and 1.644 /im transitions 
toward HH54E+K, whilst for HH7 we adopted the estimate E{J — K) — 0.7 obtained by 
Gredel (1996) for the nearby source HH8. 

In Figure 11, based upon the hne intensities tabulated in Table 2, we present rotational 
diagrams for each of the four circular regions for which we computed an average spectrum 
(§3.2). Two features are common to all observed positions. First, as discussed further in 
§4.4 below, the rotational diagrams all show the "zigzag" patterns characteristic of non- 
equilibrium ortho-to-para ratios, in which a fit to the para- (even- J) states lies above that 
for the ortho- (odd- J states). Second, the rotational diagrams all exhibit a degree of curva- 
ture, with the slope of the rotational diagram decreasing with increasing J. This curvature 
indicates that each beam is sampling material at a range of temperatures. Furthermore, 
the absence of any turnover in the curvature at high J suggests that the rotational level 
populations are close to local thermodynamic equilibrium (LTE) even for J as high as 9. 

For each circular region in Table 2, we determined an excitation temperature, T^y, from 
the ratio of the J — X and J — Y column densities and defined according to 

rp _ Ey — Ex 

k\Ti[NxgYlNYgx\ 

where Ej and gj are the energy and degeneracy of state J. Because the ortho-para ratio 
shows departures from equilibrium, the quantity Txy is most useful when X and Y are either 
both odd or both even. The excitation temperatures T42, T53, Tg4, T75, Tgg and T97 are given 
in Table 3; the inequalities T42 < Tq^ < Tge and T53 < T75 < T97 observed in every case are 
a simple manifestation of the curvature apparent in the rotational diagrams. 

We have obtained two-component fits to the rotational diagrams, in which a least squares 
fit was obtained for a combination of warm and hot gas components, assumed to be in LTE 
except with non-LTE ortho-to-para ratios, with temperatures T^j and (typically found 
to be ~ 400 K and ~ 1000 K, respectively). The solid hnes in Figure 11 show the two- 
component fits to the rotational diagrams, while the points indicate the measured values. 
The two-component fits provide an excellent fit to the data, although the solution is clearly 
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non-unique. Indeed, there is no specific evidence for a bimodal temperature distribution, 
and the data could be fit equally well by models that invoke a continuous range of gas 
temperatures. The assumption of LTE is clearly valid for the lower rotational states, but 
the data would allow modest departures from LTE for the highest states that we have 
observed. Indeed, calculations by Le Bourlot et al. (1999) indicate that the critical density 
for J = 9 (at which the probabilities of coUisional and radiative de-excitation are equal) is 
~ 4 X lO^cm"^, a factor of several larger than the densities inferred from observations of 
CO in some of these sources (see discussion below). Departures from LTE, if present, would 
be an increasing function of the rotational quantum number, J, and would therefore require 
higher gas temperatures than the values inferred under the assumption of LTE. 

The best-fit temperatures, computed under the assumption of LTE, together with the 
beam-averaged column densities {Nyj and N^) and the H2 ortho-to-para ratios (OPR^y and 
OPR^), are hsted in Table 3. Table 3 also shows the sum of the A^(4) though A'"(9) column 
densities, A^4_9, all of which are derived directly from observations with the Short-Low mod- 
ule; not surprising, A^4-9 falls a factor of several below (A^^ + Nh) because it only includes 
states with J > 4. The estimated gas temperatures, and Th, are characteristic of low- 
velocity nondissociative shocks (which are of "C"-type in the designation of Draine (1980); 
see also Appendix B). By means of a simple analytic treatment of energy conservation in 
such shocks (Appendix B), constrained by the results of the detailed shock models of Kauf- 
man & Neufeld (1996; hereafter KN96), we find the characteristic gas temperature within 
such shocks to be given by 

7; = 375r°-^^lf K 

(eqn. B8), where lOi'sekms"^ is the shock velocity and 6(no/cm~^)^/^ yuG is the assumed 
preshock magnetic field for a shock propagating in material of preshock H2 density no- Given 
our standard assumption about the magnetic field, b — 1, the warm and hot gas components 
typically correspond to shock velocities of ~ 10 and ~ 20 kms~^ respectively. The H2 
rotational diagrams predicted by KN96 and other theoretical shock models (e.g. Wilgenbus 
ct al. 2000) all exhibit less curvature than our observed rotational diagrams, implying the 
presence of a range of shock velocities within each of the circular apertures demarked in 
Figure 7. 

The preshock H2 density, Uq, is less well constrained by the H2 line fluxes. The best 
estimates for the density of the warm shocked gas have been obtained by observations of 
high-lying rotational lines of CO, for which non-LTE line ratios provide a valuable probe. 
Based upon (9 observations of such transitions toward HH54, Liseau et al. (1996) estimated 
the H2 density in the warm shocked gas as 1 — 4 x 10^ cm""^. Since most of the high- 
J CO emission occurs before the compression of the neutral species becomes large, the 
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corresponding preshock density is only a factor ~ 1.5 smaller^ (Appendix B, eqn. B6): 
77-0 = 6 — 25 X 10^ cm~^. Similar values were obtained for HH7-11, based upon measurements 
of CO line ratios (Molinari et al. 2000), although measurements of the CO/H2 line ratios 
argued for smaller values (no ~ 10^ cm~^) toward HH7. 

The observed column density of shocked H2 is predicted to be 

Ns = 4x10^° bn°fvj-''^s, 

(Appendix B, eqn. B7), where $5 is a geometric factor defined to be the ratio of the projected 
beam size to the surface area of the shock. In principle, the quantity $s can be either larger or 
smaller than unity. For a planar shock with a beam covering factor of unity that is viewed at 
angle 6 to the shock normal, the observed vahie of <l>s is sec 6 > 1. On the other hand, $s can 
be smaller than unity when the beam is incompletely covered by shocked material. For our 
two component fits to the circular apertures presented in Table 3, the values of $s required 
for the warm and hot components are typically 0.2 n^^'^ and 0.05 n^^'^, respectively, implying 
that the shocked material has a covering factor considerably smaller than unity and possesses 
substructure that is not resolved by these observations. In HH7, the subarcsecond resolution 
maps of near-IR H2 vibrational emissions presented by Norciga-Crespo et al. (2002) clearly 
indicate that $s ^ 1 for observations with the IRS instrument. 

The spectral fine maps obtained for HH54 and HH7 also permit us to obtain a ro- 
tational diagram at each spatial position, and thereby to derive detailed maps of the gas 
temperature, ortho-to-para ratio, and column density of the shocked gas. While the aver- 
age spectra presented in Figure 8 are of sufficient signal-to-noise ratio (SNR) to permit the 
two-component fits described above, the lower SNR obtained in a single spectrum does not 
warrant a two-component fit to the rotational diagram. Accordingly, in constructing maps 
of the gas temperature, we obtained a least-squares fit to the A^(4) through A^(7) column 
densities (all measured with the SL module) for a single gas component with a temperature 
and H2 ortho-to-para ratio, maps of which are presented in Figure 12 (middle and lower 
panels respectively). Here, the white contours show the H2 S(3) line intensity, and regions 
of low SNR have been masked and appear in black. The upper panel in Figure 12 shows 
the sum of the A^(4) though A^(9) column densities, A^4-9. Figure 12 shows that the fitted 
gas temperature is remarkably constant over the mapped region, lying generally within the 
range 900 ± 200 K, and indicating that a similar admixture of shock velocities is present 
along each line-of-sight. The ortho-to-para ratio, however, shows substantial variations - as 
noted previously by Wilgenbus et al. (2001) - and is discussed further in §4.4 



^Liseau et al. significantly overestimated the compression of the neutral material, thereby obtaining a 
significant underestimate of the preshock density required to match the CO line ratios. 
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4.3. Upper limits on the H2O abundance 

The increase in temperature and density behind a shock front can alter significantly 
the composition of the gas in its wake. Specifically, chemical reactions not favored at tem- 
peratures characteristic of the preshock quiescent material (i.e., T ^ 50 K) can proceed 
rapidly in the warm post-shocked gas. One such example are the reactions that lead to 
the production of gas-phase H2O. In the cold preshock gas, cosmic-ray driven ion-neutral 
reactions are capable of converting O and H into water yielding a peak water abundance 
relative to H2, a;(H20) = n(H20)/n(H2), of ~ 10~^ in approximately 10^ years (Bergin et 
al. 1998). When the gas temperature exceeds ~300 K, a series of neutral-neutral reactions, 
H2 + O ^ OH + H and H2 + OH H2O + H, dominate and chemical models (e.g., Elitzur 
& de Jong 1978; Neufeld et al. 1995) predict that most of the gas-phase oxygen not bound 
as CO will be processed into water vapor. These reactions are rapid in warm gas, producing 
x(H20) ~2xl0~^ in less than a few hundred years (Bergin et al. 1998), but are negligibly 
slow at low temperatures because they possess significant activation barriers. 

Observations of both HH54 and HH7 have been conducted using ESA's Infrared Space 
Observatory (ISO). ISO^s Long Wavelength Spectrometer (LWS) was used in its grating 
mode to acquire low resolution (A/AA ~200) spectra between 43 and 197 //m within an 
80" diameter beam. Toward HH54, Liseau et al. (1996) detected a number of low-lying 
(i.e., E/k ^ 300 K above the ground state) H2O and OH fines, along with CO emission 
lines from the rotational upper level J = 14 to 19. These data are best fit by a model 
that assumes collisional excitation behind a 10-15 km s~^ nondissociative C-type shock and 
x(H20) =10-^ 

Toward HH7, Molinari et al. (2000) find evidence for the presence of both fast dissocia- 
tive J-type and slower C-type shocks. In particular, Molinari et al. find that a 15-20 km s~^ 
C-type shock accounts well for the observed strengths of the CO J = 14-13 to J = 17- 
16 transitions plus the H2 (0-0) S(0)-S(5) transitions, which were measured with the ISO 
Short Wavelength Spectrometer. Assuming that the ISO-detcctcd H2O 179.5 fim 2i2 — Iqi 
and 174.6 /xm 3o3— 2i2 lines arise predominately from behind the same C-type shock, Moli- 
nari et al. derive a;(H20) <6xl0~^. Mofinari et al. suggest that the lower-than-expected 
post-shock H2O abundance is due to the rapid condensation of water onto dust grains in the 
post-shocked gas. A similar effect is observed toward IC443, in which a low water abundance 
(2xl0~^ < a;(H20) <3xl0~^) is inferred^ in the presence of both C- and J-type shocks 
(Snell et al. 2005). Snell et al. ascribe the low water abundance to a combination of two 



•^Snell et al. expressed their findings in terms of the ratio n(o — H20)/n(C0). Here we make the assump- 
tion that n(C0)/n(H2) = IQ-*. 
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effects. First, more tfian fialf of the preshock oxygen not in CO remains trapped as water-ice 
on grain surfaces despite the passage of a weak C-type shock, suppressing the production 
of gas-phase H2O in the post-shocked gas. Second, ultraviolet radiation generated by the 
nearby fast J- type shock photodissociates much of the gas- phase water produced behind the 
C-type shock. 

Eight H2O transitions with upper state energies Eu/k < 1000 K lie within the wave- 
length range accessible to Spitzer/IR'^. Prior to the measurements reported here, none of 
these transitions had been obscrvationally investigated toward either HH54 or HH7. To as- 
sess the extent to which the upper limits to the H2O line strengths set here constrain the 
water abundance, the equilibrium level populations of all ortho and para rotational levels of 
the ground H2^^0 vibrational state with energies E/k up to 7700 K have been calculated us- 
ing an escape probabihty method described by Neufeld & Melnick (1991). Only two updates 
have been made in the calculations undertaken here: (1) for coUisionally induced transitions 
among the lowest 45 levels of ortho- and para-water, the rate coefficients calculated by Green, 
Maluendes, & McLean (1993) for He-H20 collisions were used (multiplied by 1.348 to ac- 
count for the different reduced mass when H2 is the collision partner), and (2) the expression 
for the photon escape probability used by Neufeld & Kaufman (1993) was adopted. 

We obtain level populations for a plane-parallel emitting region in which the limit of 
large velocity gradient applies. A two-temperature component model is used, consistent with 
the results of our H2 measurements (which arc summarized in Table 3). In addition, preshock 
densities of 10*^ and 10^ cm^'^ arc considered. The H2 column densities and velocity gradients 
are assumed to be those consistent with the simple shock model presented in Appendix B, 
i.e.: 



N{H2) = 7.9 X 10' 



,20 





0.5 




-0.555 


105 




1000 





cm 
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dv 
dz 



1.6 X 10^ 





0.5 




1.30 


105 




1000 





km s ^ pc ^, 



where n{H2) is the preshock gas density and Tgas is the gas temperature. Finally, we assume 
that the fraction of the beam filled by the warm and hot H2O is given by the ratio of the 
measured values of A'"(H2) (given in Table 3) to those derived above. For these conditions, 
the strongest water lines within the Spitzer/IK^ wavelength range are predicted to be the 
29.837 jim 725-616 and 30.899 jivn 634-605 ortho-H20 transitions. Table 4 hsts the Spitzer- 
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measured 3a upper limits to the line intensities for each transition along with the derived 
upper limits to the gas-phase water abundance. 

Several things are apparent. First, the upper limits to the H2O abundance given in 
Table 4 are higher than, but nonetheless consistent with the water abundances derived 
previously from the ISO data. In particular, for a preshock density of 10^ cm~^ - favored by 
the H2 data toward HH54 (Liseau et al. 1996) - the inferred abundance limits arc less than 
would be expected from the full conversion of elemental oxygen not locked in CO into H2O. 
Unfortunately, the present 29.837 and 30.899 /xm line flux limits are insufficient to allow us 
to comment further upon the previous HH54 models. Second, the H2O 174 and 179 /xm 
line fluxes measured by Liseau et al. provide an additional constraint on a;(H20) within the 
warm and hot components discussed here. In particular, the /jSO-measured flux in the H2O 
174 jMn. line, ~3xl0~^° W cm~^, requires that x(H20) < 2xl0~^ in the warm and hot H2 
components or these components will produce a 174 //m line flux in excess of that measured. 
This result is in agreement with the water abundance determined by Liseau et al. of 10~^ 
and suggests that the size of the H2O 174 and 179 /im emitting region is larger than that 
inferred for the S'^^ifoer-measured H2 emission (or that assumed for the mid-infrared H2O 
emission). This is plausible since the 174 and 179 fim lines probe cooler gas than either the 
H2 or mid-infrared H2O lines. 

Toward HH7, Molinari et al. estimate the preshock density to be ~ IC^ cm^'^. For this 
density, rather high H2O abundances are required to reproduce the 29.837 and 30.899 ^m. 
line intensity limits. However, the strength of the H2O 179 ^m. line observed by Molinari 
et al. does set a more stringent upper limit of a;(H20) <6xl0^"^ within the warm and hot 
components measured here. Unfortunately, this limit remains significantly greater than the 
water abundance inferred from the ISO measurements and, thus, does not permit a test of 
their model. 



4.4. The H2 ortho-to-pEira ratio 

The H2 ortho-to-para ratios (Table 3) derived for the two-component fits, OPR^ and 
OPR^, are uniformly less than the local thermodynamic ratio, OPRlte ~ 3 at the tempera- 
tures of relevance. This non-equilibrium behavior, first observed in the atmospheres of giant 
planets (e.g. Conrath & Gierasch 1983), was noted in ISO observations of HH54 (NMH98; 
Wilgenbus et al. 2001) and subsequently observed in several other interstellar regions (e.g. 
Fuente et al. 1999, Rodriguez-Fernandez et al. 2000). As discussed in NMH98, the H2 ortho- 
to-para ratio equilibrates very slowly, exhibiting "memory effects" in which the OPR can 
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reflect the temperature of a previous epoch in its thermal history^. Thus NMH98 argued 
that the OPR of 1.0 ± 0.4 observed toward HH54 indicated (1) that the warm (T ~ 650 K) 
shocked H2 that they observed in HH54 must have previously acquired an OPR character- 
istic of equilibrium at a temperature < 90 K and (2) that the observed gas has not been 
warm long enough for a high-temperature OPR to have been established. 

The equilibration of ortho- and para-H2 has been considered in the detailed shock models 
of Timmermann (1998) and Wilgenbus et al. (2000). Both studies concluded that reactive 
collisions with atomic hydrogen are the dominant para-to-ortho conversion process in non- 
dissociative molecular shocks: 

para-H2 + H ^ H + ortho-H2 (Rl) 

Because reaction (Rl) has a substantial activitation energy barrier (AEA/k ~ 3900 K), the 
efficiency of para-to-ortho conversion is a strong function of temperature (and, consequently, 
shock velocity). This strong temperature dependence is entirely consistent with the results 
in Table 3, which show that for each region the hot gas component has a much higher OPR 
than the warm gas component. This behavior is also apparent in the H2 rotational diagrams 
(Figure 11), in which the departures from LTE (i.e. the degree of zigzag) are plainly larger 
for small J than for large J. 

In Figure 13, the ortho-to para ratios, OPR^y and OPR^, for the two-component fits in 
Table 3, are plotted as a function of the fitted gas temperatures, and Th (open triangles). 
We have also obtained two-component fits to the gas temperature and OPR elsewhere in 
the mapped regions. Averaging the observed spectra in a set of 5" x 5" subregions to attain 
adequate SNR to yield meaningful two-component fits, we determined OPR^,,, OPR/^, 
and T/j estimates for each such subregion. Green and red squares in Figure 13 indicate 
respectively the (T„,,OPR„,) and {Th, OPKh) values obtained in subregions within the HH54 
and HH7 maps. Green and red crosses apply to the discontiguous regions northwest of 
HH7 (Figure 12; these are close to HH8). The orange curve denotes the ortho-to-para ratio 
expected in LTE. 



^Such effects are a well-known problem in the industrial production and storage of liquid hydrogen (LH2). 
The rapid refrigeration of molecular hydrogen leads to the production of LH2 with an ortho-to-para ratio 

3, reflecting the "frozen-in" equilibrium value at room temperature, and greatly exceeding the value of 
0.02 obtained at equilibrium at 21 K, the boiling point of LH2. Equilibration occurs over a timescale of 6.5 
days, resulting in the eventual release of a total heat content that exceeds the latent heat of vaporization of 
LH2: without continued refrigeration, the LH2 simply boils away! The solution to this problem involves the 
introduction of a paramagnetic catalyst to facilitate ortho-to-para conversion during the original liquification 
process and thereby prevent the room temperature OPR from being frozen-in. 
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The green and red symbols in Figure 13 confirm that the ortho-to-para ratio tends to 
increase with gas temperature. The black curves show the expected behavior for a simple 
model in which all the observed gas started with the same initial ortho-to-para ratio - 
OPRq = 0.4 (corresponding to equilibrium at temperature 55 K) for HH54 and OPRq = 0.25 
(corresponding to equilibrium at temperature 48 K) for HH7 - and has been at the currently- 
observed temperature T for the same time period, r. In that case, the current OPR is given 



where n(H) is the atomic hydrogen density and kpo is the rate coefficient for para-to-ortho 
conversion, estimated as 8 x 10~^^ exp(— 3900/T) cm^ s~^ by Schofield et al. (1967), based 
upon the laboratory experiments of Schulz & Lc Roy (1965). The black curve applies to the 
case n2(H)r = 120 yr, where n2(H) = n(H)/10^ cm~^. If, following NMH98, we assume that 
H atoms are produced primarily by reaction of O with H2 to form OH, followed by reactions 
of OH with H2 to form water, we expect an n(H)/n(H2) ratio ~ 10~^, corresponding to 
77,2 (H) ~ 1 for a preshock H2 density Uq — lO^cm"^. The steep rise in the black curve 
at T ~ 1000 K is a consequence of the strong temperature-dependence of kpQ. To within 
the likely errors, all the points lie between the black and orange curves in Figure 13. This 
behavior then implies that all the shocked gas had an initial OPR between 0.4 and 3 and has 
been warm for at least 150 n2(H)~^ yr. The foregoing analysis is premised on the assumption 
that the rotational temperature inferred for the hot gas component is equal to the kinetic 
temperature. As discussed in §4.2, above, subthcrmal excitation of the S(7) transition is 
allowed by the data and would mean that the actual gas temperature is higher than the 
rotational temperature for the hot component, T^. Such effects, in turn, would decrease the 
required timescale for ortho-para conversion. 

For a postshock atomic H abundance ~ 10^^, the timescale of 150n2(H)^^ yr inferred 
above is broadly consistent with the time period for which gas temperature is elevated upon 
passing through a nondissociative shock. For uq — lO^cm"^ and Vg — 15kms~^, the simple 
treatment of Appendix B implies a shock column density, Ns ~ 10^^6 Aj„cm~^, corresponding 
to a timescale Ns/{noVs) — 200 6Aj„ yr. However, the OPR values plotted in Figure 13 are 
entirely inconsistent with a shock model with a large H abundance, as has been invoked by 
Smith et al. (2003), who proposed an n(H)/n(H2) ratio of 3 for the bow shock in HH7. Such 
a large H abundance would lead to very rapid para-to-ortho conversion, producing OPR 
values significantly in excess of those observed here. 

We may also compare the results plotted in Figure 13 with the predictions of detailed 
theoretical models for para-to-ortho conversion in non-dissociative molecular shocks. The 
cyan loci in Figure 13 represent the predictions of Wilgenbus et al. (2000) for a shock 



by 





1 + OPR(t) 1 + OPR, 
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propagatiiig in gas with an initial OPR of 0.01. The loci correspond to preshock H2 densities 
of 10^, 10^, 10^ and 10^ cm~^, respectively from left to right, and the asterisks apply to shock 
velocities of 30 and 40 kms~^ (10^cm~^ case, from lower left to upper right); 20, 30, and 
40 kms~^ (10^ and lO^cm"^ cases, from lower left to upper right); or 10, 15, 20, 25, and 
30 kms~^ (10^cm~^ case, again from lower left to upper right). As in the case of the toy 
model (black curve), the curves for 10^ or lO^cm"^ are a reasonable approximation to the 
lower envelope of the observed data points. The absence of sight-lines with high temperature 
and low OPR is consistent with shock models in which para-to-ortho conversion is driven by 
reactive collisions with atomic hydrogen. On the other hand, the presence of sight-lines with 
low temperature and high OPR suggests that some of the shocked material has a high initial 
OPR. Such material is clearly present for the sight-lines represented by the crosses in Figure 
13. These regions lie behind the "working surface" of the jet in HH7-11; they may, perhaps, 
have been heated by relatively fast shocks during a previous epoch of shock activity, thereby 
acquiring a larger OPR, and are currently being reheated by a slower shock. The bow shock 
in HH7, by contrast, appears to be propagating into fresh material in which the OPR is 
^ 0.25, corresponding to equilibrium at temperature ~ 50 K. 



5. Summciry 

We have carried out spectroscopic observations toward the Herbig-Haro objects HH7- 
11 and HH54 over the 5.2 — 37/im region using the Infrared Spectrograph of the Spitzer 
Space Telescope, and have thereby mapped the S(0) - S(7) pure rotational lines of molecular 
hydrogen, together with emissions in fine structure transitions of Ne"*", Si"*", S, and Fe"*". 

The primary results of our study are: 

1. The H2 rotational emissions indicate the presence of warm gas with a mixture of 
temperatures in the range 400 — 1200 K, consistent with the expected temperature behind 
nondissociative shocks of velocity ~ 10 — 20kms~^. 

2. The fine structure emissions are spatially offset from the H2 emissions, lying closer 
to the source of the outflow that drives the shocks in these sources. They originate in faster 
shocks - of velocity ~ 35 — 90kms~^ - that are dissociative and ionizing. Such shocks likely 
result from the deceleration of the outflow, and are offset from a slower non-dissociative bow 
shock that is driven into the ambient molecular medium. 

3. Maps of the H2 line ratios reveal little spatial variation in the typical admixture 
of gas temperatures in the mapped regions, but show that the H2 ortho-to-para ratio is 
quite variable, typically falling substantially below the equilibrium value of 3 attained at the 
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measured gas temperatures. 

4. The non-equilibrium ortho-to-para ratios are characteristic of temperatures as low as 
~ 50 K, and are a remnant of an earlier epoch, before the gas temperature was elevated by 
the passage of a shock. Correlations between the gas temperature and H2 ortho-to-para ratio 
show that ortho-to-para ratios < 0.8 are attained only at gas temperatures below ~ 900 K; 

this behavior is consistent with theoretical models in which the conversion of para- to ortho- 
H2 behind the shock is driven by reactive collisions with atomic hydrogen, a process which 
possesses a substantial activation energy barrier {E^/k ~ 4000 K) and is therefore very 
inefficient at low temperature. 

5. The lowest observed ortho-to-para ratios of only ~ 0.25 suggest that the shocks in 
HH54 and HH7 are propagating into cold clouds of temperature ^ 50 K in which the H2 
ortho-to-para ratio is close to equilibrium. 

6. Upper limits on the 29.837 /im 725-616 and 30.899 /xm 634-505 transitions of ortho-H20 
place upper limits on the H2O abundance that are higher than - and therefore consistent with 
- abundances derived previously from ISO observations of transitions at longer wavelength. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under a NASA 
contract. D.A.N, and P.S. gratefully acknowledge the support of grant NAG5-13114 from 
NASA's Long Term Space Astrophysics (LTSA) Research Program. 
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Appendix A: data reduction procedures 

Our data reduction procedure involved the following steps: 

1) Removal of bad pixels in the "Basic Calibrated Data (BCD)" files that are generated 
by the Spitzer Science Center for the Short-High and Long-High modules. Bad pixels were 
defined as those flagged in either the SSC "bmask" flle or in a "grand rogue mask" generated 
at University of Rochester from an analysis of dark fields acquired at the start of each IRS 
campaign. As a rogue might not appear in every campaign dark field, we assumed that any 
pixel identified as a rogue in ANY campaign would be corrected for all campaigns; thus, the 
grand rogue mask contains 25% more rogue pixels than a single campaign rogue mask. 
They were removed with the "imclean" routine, which interpolates between nearby good 
pixels. 

2) Eoctraction of the individual spectra at each position along the slit. These spectral 
extractions were accomplished using an augmentation to the SMART software package that 
has been developed at University of Rochester. This augmentation generates a spectrum for 
each resolution element along the slit, resulting in a series of spectra at defined positions 
in R.A. and declination. The spectra were trimmed to eliminate overlap between different 
orders of the IRS. 

3) Regridding of the spectra onto a grid that is regular in R.A. and declination. Because 
there are positions that arc multiply observed with more than one slit position, this regridding 
sometimes involves the averaging of spectra obtained at nearly identical positions. 

4) Application of the slit loss correction function (SLCF). Because the SSC- provided 
data arc flux-calibrated using point source calibrators, a downwards correction must be 
provided to obtain intensities for extended sources, the correction factor (i.e. the "SLCF") 
being the fraction of the point source flux that falls within the slit width. We used the 
SLCFs kindly provided by J. D. Smith from an analysis of the instrumental point spread 
function. 

5) Determination of line and continuum fluxes at each spatial position within the mapped 
region. Here, we used a Levenberg-Marquardt (L-M) fltting routine simultaneously to flt a 
Gaussian line and second order polynomial baseline at each position. Line fits were obtained 
for a series of spectral fines identified in the co-averaged spectra for each field. In fitting 
the co-averaged spectra, which are of high signal-to-noise ratio (SNR), the line centroids 
and widths were treated as free parameters to be optimized in the L-M fit. The resultant 
best-fit centroids and widths were then used as fixed parameters in the L-M fits to the lower 
SNR spectra obtained at individual spatial positions; for these, only the line intensities 
and continuum parameters were optimized. This procedure is justified by the fact that the 
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expected Doppler shifts associated with gas motions in the source are very small compared 
to the spectral resolution, A/ A A being only ~ 60 and ~ 600 respectively for the low (SL) 
and high (SH and LH) resolution modules. 

The data reduction procedure was accomphshed using the SMART software package 
and a set of associated IDL procedures that we have written. The contour maps presented 
in Figures 1-7 were also created using IDL procedures. 

To verify our data reduction procedure, we have conducted several tests. These include: 

1) Continuum mapping of SVS13. Figure 14 shows contour maps obtained for several 
continuum wavelengths for HH7-11 Field 1 (see Table 1), which contains the strong contin- 
uum source SVS13 at its center (i.e. at position (0,0) in Figure 14). The intensity maps 
shown in Figure 14 use logarithmic contours spaced by a factor 2^/^, with the second-highest 
contour corresponding to an intensity one-half the maximum. To the extent that SVS13 
can be regarded as an unresolved point source at these wavelengths, the maps in Figure 
14 represent the beam profile resulting at the end of the entire data reduction procedure. 
These maps allow the astrometric accuracy and the half-power-beam width (HPBW) to be 
estimated for the final data product. Table 5 shows the results of Gaussian fits to the spatial 
profiles, and indicates that the astrometric accuracy is better than ~ 2" in the final data 
product. The map at 5.8 //m clearly exhibits a sixfold rotational symmetry at intensities lev- 
els less than ~ 3% of the maximum, a behavior that is also evident in Infrared Array Camera 
(IRAC) observations (Hora et al. 2004) and reflects the diffraction pattern produced by the 
telescope. 

2) Spectral eoctraction of SVS13 and comparison with Infrared Space Observatory (ISO) 
observations. As a check of the flux calibration, we have extracted spectra of the strong 
continuum source SVS13, obtaining a co-average from the entire map of HH7-11 Field 1. For 
comparison, we have obtained an archival spectrum of SVS13 from the ISO Short Wavelength 
Spectrometer (SWS), co-averaging multiple scans and removing bad pixels with the use of 
the ISAP software package. Figure 15 compares the spectrum obtained with Spitzer (heavy 
black line) with that obtained with ISO (red, blue, green and magenta lines, different colors 
being used to show the different SWS orders). Overall, the agreement is very good, with 
the absolute fluxes typically agreeing to within ~ 15%. In the spectral region longward of 
27.5 //m, the fluxes obtained with ISO/SMVS exceed those measured by Spitzer by as much 
as 25%. However, a significant order mismatch at 27.5 //m (between detectors 11 and 13) 
in the ISO spectrum suggests that the /.SO-determined flux is an overestimate longward of 
27.5 /im. 



Appendix B: simple analytic treatment of nondissociative C-type shocks 
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Here we present a simple analytic treatment of nondissociative C-type shocks, which 
is limited in its applicability to steady-state shocks with plane parallel geometry. Detailed 
models of planar nondissociative C-type shocks (e.g. Kaufman & Ncufcld 1996, hereafter 
KN96; Wilgenbus et al. 2000), predict line strengths for the H2 pure rotational transitions 
that yield rotational diagrams with remarkably little curvature. This behavior suggests 
that the warm gas behind such shocks can be approximated as an isothermal, isobaric slab. 
In this appendix, we show that such a description yields predictions for the emission line 
spectrum of the shocked gas that are in good agreement with those of detailed models, and 
we derive analytic expressions for the appropriate slab thickness, density, temperature and 
velocity gradient as a function of the preshock H2 density, tiq = Uq^ 10^ cm shock velocity, 
Vs — 10 1)56 km s"^, and preshock magnetic field, Bq — 6 (no/cm ~^)^/^ //G. 

Nondissociative C-type shocks have been discussed extensively in the literature (e.g. 
Draine & McKee 1993, and references therein) and will not be reviewed here. They were 
originally defined by Draine (1980), who distinguished between J- (or "jump")-type shocks 
in which the fluid velocity changes suddenly and the neutral and ionized species share a 
common velocity, and C- or ( "continuous" ) -type shocks in which the fluid velocity varies 
slowly. The basic physics of C-type shocks had been discussed previously by Mullan (1971), 
who recognized ion-neutral drift to be an essential feature of magnetic shocks propagating 
in weakly-ionized media. Ions entering the shock, being strongly coupled to the magnetic 
field, are decelerated fairly rapidly and quickly reach their maximum compression. Neutral 
species entering the shock, by contrast, are uncoupled to the magnetic field and are slowly 
decelerated in collisions with the ions through which they drift. Thus the maximum heating 
rate within the shocked material occurs where the compression of the ions is large and that of 
the neutral species is small. Accordingly, we make the following simplifying approximations: 

1) The neutral fiuid is at a constant compression, C„ = n(H2)/no, which is slightly 
greater than unity; the exact value is an adjustable parameter to be determined below by a 
comparison with the fine fiux predictions of the KN96 shock models. 

2) The ionized fiuid is at a constant compression equal to the maximum that is attained 
behind the shock 

Ci = ni/mo = 2 M'f = 7.7vsG/b, (Bl) 

where Ma — Vse/iO.lSAb) is the Alfven Mach number and rijo is the preshock ion density. 
Additional ionization within the shock ( "self-ionization" ) is neglected. 

3) The thickness of emitting slab is given by 

t = kiLinino, Vs)/Ci, {B2) 
where Lj„ (n[H2], I'd) is the ion-neutral couphng length (KN96) for H2 density n(H2) and ion- 
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neutral drift velocity Vd, and ki is a dimensionless constant of order unity, the exact value 
of which is an adjustable parameter to be determined below by a comparison with the line 
flux predictions of the KN96 shock models. 

4) The slab is at constant temperature, T, determined such that the radiative output is 
equal to the mechanical luminosity, implying 

nlCltA{T, noCn) = lfinov',{l - [2Ma]-^) ~ IfiUovl (S3) 

where A(T, n[H2]) is the cooling rate coefficient for temperature, T — lO^TaK, and H2 
density, n{R2) = lO^nscm"^. 

5) The velocity gradient is dv^/dz = k^Vs/t, where k'2 is a dimensionless constant of 
order unity, the exact value of which is an adjustable parameter to be determined below by 
a comparison with the line flux predictions of the KN96 shock models. 

6) We adopt as our standard assumptions an ion-neutral coupling length 

L^„ = 3.5xlOX-°-^^°fcm, {BA) 

based upon a fit to results of KN96 (obtained for the case where polycyclic aromatic hydro- 
carbons (PAHs) are assumed to be present), and a cooling rate coefficient 

A°(r, n[H2]) = 1.4 X lO-'^ng-^-^Tl'^ ergcm^ s"^ (S5) 

based upon a flt to the results of Neufeld et al. (1995). The superscripts on Lf^ and A° are 
used to denote the "standard values" adopted for our treatment. 

Assumptions (1) - (6) above define the required properties of the isothermal, isobaric 
slab, once the three parameters Cm ki, and /c2 have been specified. In determining the best- 
fit values for these parameters, we have compared the emission predicted for the slab model 
with the detailed predictions of KN96. Here we assumed fixed CO and H2O abundances, 
corresponding to the assumption that CO and H2O account for all the O and C nuclei in the 
gas-phase. Because the O and C abundances assumed by KN96 were somewhat larger than 
those adopted by Neufeld ct al. (1995) in computing the cooling rate coefficient, A°(T, n[H2]), 
we adopted A(T, n[H2]) = 1.5 A°(T, Ti[H2]) for the purposes of this comparison. The optimum 
fit to the KN96- predicted emission line spectrum is obtained for C„ = 1.5, /ci = 2, and k2 = 2. 

Combining the expressions given above, we find that the slab H2 density is fit best by 

n(H2) = 1.5 X 10^ no5 cm'^, (S6) 

the H2 column density by 

A^(H2) = 1.37 X 10^1 n',iv:,'-'' b {hjLl) cm'', {B7) 
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the gas temperature by 

T = 375 vlf (L,„/L°J-o-36 (A/AO)-o-36 k (^8) 
and the velocity gradient by 

dv/dz = 4500 n'^oi b'^ {1^^/ L'^J'^ kms"^ pc"^ (59) 

In Figure 16, we compare the H2 S(5)/S(3) hne ratio and the total H2, H2O and CO 
surface luminosities derived for the simple slab model (red) with those obtained in the 
detailed calculations of KN96 (blue). In general, the agreement is remarkably good. Almost 
invariably, either (a) the shock velocity corresponding to a given flux or flux ratio agrees to 
within ~ 5kms~^; and/or (b) the predicted line flux or flux ratio for a given shock velocity 
agrees to within a factor 2. Often, the agreement is considerably better. The main exception 
is for water emissions at low shock velocity (vg < 15 kms""*^), where the detailed KN96 models 
predict considerably less H2O than the slab models. This behavior results because - for the 
slab models - we assumed that water vapor accounts for all gas-phase oxygen nuclei not 
locked in CO, whereas the treatment of oxygen chemistry in KN96 predicts considerably 
lower water abundances at low shock velocity. 
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Table 1. Details of the observations 



Source 


Date 


Module 


(Jbservmg 


Map center 


Map size 








time (s) 


R.A. (J2000) 


Dec. (J2000) 


(arcsec) 


HH7-11 Field 1 


2004 Feb 15 


SL 




3n zym 3. i Is 


+31a lb 01. y 


58 X 57 


HH7-11 Field 1 


2004 Feb 15 


SH 




3h 29m 3.66s 


+31d 16' 03.6" 


58 X 76 


HH7-11 Field 1 


2004 Feb 15 


LH 


10243^ 


31i 29m 3.66s 


+31d 16' 03.6" 


61 X 77 


HH7-11 Field 2 


2004 Feb 16 


SL 




31i 29m 8.61s 


+31d 15' 26.6" 


58 X 57 


HH7-11 Field 2 


2004 Feb 16 


SH 




31i 29m 8.56s 


+31d 15' 28.3" 


58 X 76 


HH7-11 Field 2 


2004 Feb 16 


LH 


10246^ 


31i 29m 8.56s 


+31d 15' 28.3" 


61 X 77 


HH54 


2004 Apr 7 


SL 


1133 


12h 55m 51.52s 


-76d 56' 19.1" 


47 X 64 


HH54 


2004 Apr 7 


SH 


3260 


12h 55m 50.87s 


-76d 56' 19.3" 


47 X 52 


HH54 


2004 Apr 7 


LH 


931 


121i 55m 50.82s 


-76d 56' 18.0" 


44 X 67 



^Total for SL, SH, and LH modules 
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Table 2. Beam-averaged line intensities 



Transition 


Wavelength 


Intensity (10 ^ 


erg cm s ^ sr ^ ) 






HH54FS 


HH54C 


HH54E+K 


HH7 


H2 8(0) J = 2-0 


28.2188 


0.52 


0.44 


0.47 


0.78 


H2 S(l) J = 3-l 


17.0348 


3.68 


3.18 


3.99 


3.07 


H2 S(2) J = 4-2 


12.2786 


14.1 


21.4 


17.2 


18.2 


H2 S(3) J = 5 - 3 


9.6649 


24.7 


20.7 


29.5 


20.6 


H2 S(4) J = 6 - 4 


8.0251 


24.5 


28.9 


22.9 


25.6 


H2 S(5) J = 7 - 5 


6.9095 


33.0 


24.4 


32.9 


29.0 


H2 S(6) J = 8 - 6 


6.1086 


17.0 


14.0 


14.3 


12.0 


H2 S{7) J = 9 - 7 


5.5112 


15.6 


12.7 


17.2 


16.0 


Nc+ P" 
^^•^ -^1/2 -^3/2 


12.8149 


1.68 


1.03 


0.82 


< 0.36 (3o-) 


S ^Pi -3 P2 


25.2490 


1.75 


0.75 


0.75 


1.73 


0;+ 2 pO 2 pO 

-^3/2 M/2 


34.8141 


6.07 


2.95 


3.43 


4.18 


Fe+ ^9/2-609/2 


5.3403 


8.6 ± 0.8" 


6.9 ±0.8 


2.6 ± 1.0 




Fe+ *F7/2-4F9/2 


17.9363 


0.88 ±0.18 


0.55 ±0.08 


< 0.5 (3(t) 


< 0.4 (3o-) 


Fe+ 4f5/2_4f,/2 


24.5186 


0.23 ± 0.04 


0.16 ±0.02 


0.19 ±0.03 


< 0.12 (3o-) 


Fe+ 6D7/2-®D9/2 


25.9882 


3.33 


1.69 


1.38 


1.69 


Fe+ «D5/2-«D7/2 


35.3491 


0.92 ±0.17 


0.45 ±0.11 


0.56 ±0.08 


0.57 ±0.06 


'^Error bars, where g 


;iven, are lu statistical errors. 


When error 


bars are not i 


given, the uncer- 



tainty is dominated by systematic errors, which we estimate to be < 25% (see Appendix A). 
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Table 3. H2 parameters 





HH54FS 


HH54C 


HH54E+K 


HH7 


Rotational state 


logio 


(Column density in cm 




J=2 


19.51 


19.44 


19.47 


19.70 


J=3 


18.94 


18.88 


18.98 


18.89 


J=4 


18.62 


18.81 


18.71 


18.77 


J=5 


18.25 


18.17 


18.33 


18.26 


J=6 


17.70 


17.77 


17.67 


17.75 


3=7 


17.40 


17.27 


17.40 


17.36 


J=8 


16.78 


16.69 


16.70 


16.65 


J=9 


16.45 


16.36 


16.49 


16.48 


Rotational temperatures 


(K) 








T42 


445 


571 


501 


428 


T53 


729 


717 


764 


779 


Te4 


716 


650 


647 


658 




919 


871 


852 


877 




986 


857 


946 


841 


Tgr 


1076 


1120 


1120 


1156 


Two-parameter fits 










T„ (K) 


382 


558 


470 


422 


logio (N^/cm-'^) 


19.89 


19.80 


19.84 


19.99 


OPR^ 


0.43 


0.41 


0.48 


0.21 


Th (K) 


1025 


1150 


1062 


1029 


logio (AT^/cm-^) 


19.22 


18.86 


19.13 


19.17 


OPRh 


1.53 


1.55 


1.98 


1.88 


logio (N-u, + Nh)/cra-'^) 


19.97 


19.84 


19.92 


20.05 


N4-9 (cm-2) 


18.83 


18.88 


18.91 


18.93 



-28- 



Table 4. Spitzer-Dehved H2O Abundance Upper Limits 







7 25 - 6 16 


(29.837 /im) 


634 - 5o5 


(30.899 /xm) 


Source 


H2 density (cm ^) 


Line intensity 


' n(H20)/n(H2) 


Line intensity 


' n(H20)/n(H2) 


HH54 FS 


10^ 


<9 


< 1.6 xlO-^ 


<9 


< 1.6 xlO-^ 




10^ 


<9 


< 1.7 xlO-^ 


<9 


< 1.5 xlO-=^ 


HH54 C 


10^ 


<6 


< 1.1 xlO-^ 


<4 


< 6.6 xlO-^ 




10^ 


<6 


< 1.2 xlO-^ 


<4 


<6.1 xlO-^ 


HH54 E+K 


10^ 


<7 


< 1.2 xlO-^ 


<9 


< 1.5 xlO"^ 




10^ 


<7 


< 1.3 xlO-^ 


<9 


< 1.4 xlO-=^ 


HH7 


10^ 


< 10 


< 1.8 xlO"^ 


< 14 


< 2.5 xlO"^ 




10^ 


< 10 


< 1.9 xlO-3 


< 14 


< 2.3 xlO-^ 



^3(7 upper limits in units of 10 ''erg cm ^s ^ sr ^ 
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Table 5. Effective beam size and centroid accuracy 



Wavelength 
/im 


Centroid offset 
(R.A., Dec.) 


Half-power beam width 
// 

(circular average) 


5.8 


( +0.5, -1.9 ) 


3.0 


9.2 


( -0.6, -2.3 ) 


3.1 


10.5 


( +0.1, +1.6 ) 


5.0 


18.5 


( -0.5, +2.0 ) 


5.3 


22.0 


( -0.2, +0.0 ) 


9.6 


30.5 


( -0.3, -0.4 ) 


10.2 
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Fig. 1. — H2 S(0) - S(3) emission line intensities observed toward HH54. Logarithmic 
contours are shown, in steps of 0.1 dex, with the lowest contours corresponding to line 
intensities of 3 x 10~^, 2 x 10~^, 3 x 10~^, and 1 x 10~^erg cm~^s~^sr~^ respectively for H2 
S(0), S(l), S(2) and S(3). The straight lines demark the regions within which each transition 
was mapped. The horizontal and vertical axes show the R.A. {Aa cos 5) and declination (AS) 
offsets in arcsec relative to a =12h 55m 53.40s, S = -76d56'20.5" (J2000). 
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Fig. 2. — H2 S(4) - S(7) emission line intensities observed toward HH54. Logarithmic 
contours are shown, in steps of 0.1 dex, with the lowest contours corresponding to a line 
intensity of 2 x 10^^ ergcm^^ s~^ sr~^. The straight lines demark the regions within which 
each transition was mapped. The horizontal and vertical axes show the R.A. (Act cos S) and 
declination (AS) offsets in arcsec relative to a =12h 55m 53.40s, S = -76d56'20.5" (J2000). 



20 10 -10 -20 -30 -40 20 10 -10 -20 -30 -40 




20 10 -10 -20 -30 -40 20 10 -10 -20 -30 -40 



Fig. 3. — Fine structure emission line intensities observed toward HH54. Logarithmic con- 
tours are shown, in steps of 0.1 dex, with the lowest contours corresponding to line intensities 
of 1 X 10~^, 1 X 10~^, 5 X 10~^, and 2 x 10~^erg cm~^s~^sr"^ respectively for the [Nell] 
12.8 /xm, [Fell] 26 /xm, [SI] 25 /xm and [Sill] 35 transitions. The straight lines demark the 
regions within which each transition was mapped. The horizontal and vertical axes show 
the R.A. (Aq;cos5) and dechnation (AS) offsets in arcsec relative to a =12h 55m 53.40s, 
5 = -76d56'20.5" (J2000). 
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H2 S(2) - SH H2 S(1) 




80 60 40 20 -20 -40 80 60 40 20 -20 -40 



Fig. 4. — H2 S(l) - S(3) emission line intensities observed toward HH7-11. Logarithmic 
contours are shown, in steps of 0.1 dex, with the lowest contours corresponding to line 
intensities of 5 x 10~^, 3 x 10~^, 3 x 10~^, and 1 x 10~^ erg cm"^ s~^ sr~^ respectively for 
H2 S(l), S(2) and S(3). The straight hnes demark the regions within which each transition 
was mapped, and the dashed circles centered at offset (0,0) encircle the regions close to 
the bright continuum source SVS13 where the line intensity estimates are unrehable. The 
upper map labeled H2 S(2) was obtained with the Short-High module and the lower map 
with the Short-Low module. The horizontal and vertical axes show the R.A. (Aq;cos(5) and 
declination (AS) offsets in arcscc relative to a =3h 29m 3.74s, 6 = +31dl6'04.0" (J2000). 
The marked positions in the upper two panels are based upon the observations of Khanzadyan 
et al. (2003). 
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H, :;(4) 




80 60 40 20 -20 -40 



H2 S(6) 




80 60 40 20 -20 -40 




-80 \ 1 

80 60 40 20 -20 -40 



H2 S(7) 




-80 \ 1 

80 60 40 20 -20 -40 



Fig. 5. — H2 S(4) - S(7) emission line intensities observed toward HH7-11. Logarithmic 
contours are shown, in steps of 0.1 dex, with the lowest contours corresponding to line 
intensities of 2 x 10"'^ ergcm"^ s~^ sr~^. The straight lines demark the regions within which 
each transition was mapped, and the dashed circles centered at offset (0,0) encircle the 
regions close to the bright continuum source SVS13 where the line intensity estimates are 
unreliable. The horizontal and vertical axes show the R.A. (Act cos 6) and declination (AS) 
offsets in arcsec relative to a =3h 29m 3.74s, S = +31dl6'04.0" (J2000). 
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Fig. 6. — Fine structure emission line intensities observed toward HH7-11. Logarithmic 
contours are shown, in steps of 0.1 dex, with the lowest contours corresponding to line 
intensities of 1 x 10~^, 5 x 10~^, and 2 x 10~^erg cm~^s~^sr~^ respectively for the [Fell] 
26 nm, [SI] 25 /im and [Sill] 35 /xm transitions. The straight lines demark the regions within 
which each transition was mapped. The horizontal and vertical axes show the R. A. {Aa cos S) 
and declination (AS) offsets in arcsec relative to a =3h 29m 3.74s, 6 = +31d 16'04.0" (J2000). 
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Fig. 7. — Comparison between the distribution of H2 and fine-structure emission lines. Green 
contours: H2 S(3) line; red contours [Nell] 12.8/im line (HH54) or [Sill] 35 /xm line (HH7) 
Logarithmic contours are shown, in steps of 0.1 dex, with the lowest contours corresponding 
to line intensities of 2 x 10~^, 2 x 10~^ and 3 x 10~^ erg cm~^ s~^ sr~^ respectively for the 
H2 S(3), [Nell] and [Sill] transitions. The dashed circles indicate the apertures for which 
average spectra are shown in Figure 8: these are denoted HH54E+K, HH54FS, HH54C (left 
panel from south to north) and HH7 (right panel). 
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Fig. 8. — Average spectra observed for 15" (HPBW) diameter circular apertures cen- 
tered at offsets (Aacos(5,A5) = (-7", +2") (HH54FS), (0", +16") (HH54C) and (4", -3") 
(HH54E+K) relative to a =12h 55m 53.40s, 6 = -76d56'20.5" (J2000), and offset 
(Aacos5,A5) = (+60", -37") (HH7) relative to a =3h 29m 3.74s, 6 = +31dl6'04.0" 
(J2000). The spectra in the left panel were obtained with the SL module, while those 
in the right panel were obtained with SH (A < 19.5 //m) and LH (A > 19.5 //m). 
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Fig. 9. — [Fell] fine-structure lines detected toward HH54. The spectra are averages for a 15" 
(HPBW) diameter circular aperture centered at offset {Aa cos 5, AS) — (—7", 2") (position 
HH54FS) relative to a =12h 55m 53.40s, S = -76d56'20.5" (J2000). Arbitrary continuum 
offsets have been introduced for clarity. 
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Fig. 10. — Grotrian diagram for [Fell], showing the five fine-structure transitions detected 
toward HH54 
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Fig. 11. — H2 rotational diagrams, computed for 15" (HPBW) diameter circular aper- 
tures centered at offsets {Aa cos 5, AS) = (-7", +2") (HH54FS), (0",+16") (HH54C) and 
(4", -3") (HH54E+K) relative to a =12h 55m 53.40s, 5 = -76d56'20.5" (J2000), and off- 
set {AacosS,AS) = (+60", -37") (HH7) relative to a =31i 29m 3.74s, 5 = +31dl6'04.0" 
(J2000). Square symbols indicate the measured values, while the lines show the best two- 
component fit to the observations. 
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Fig. 12. — Maps of the H2 column density (top panels), H2 rotational temperature (middle) 
and H2 ortho-to-para ratio (bottom panels) derived from a least-squares fit to the J = 4, 5, 
6, and 7 column densities. Superposed white contours (middle panels) show the H2 S(3) line 
intensity (contours at 3.0, 4.5, 6.0, 7.5 xlO~^erg cm~^ s^^ sr^^ for HH54 and 2.5, 3.5, 4.5, 
5.5 xlO~^erg cm^^ s^^ sr~^ for HH7.) Regions where the signal-to-noise ratio is inadequate 
to yield a fit to the H2 rotational diagram appear in black. 
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Fig. 13. — Correlation between ortho-to-para ratio and rotational temperature, obtained 
from two-component fits to the J = 2, 3, 4, 5, 6, 7, 8 and 9 column densities. Triangles are 
from a two-component fit to the rotational diagrams shown in Figure 11. Other data points 
apply to 5 X 5 arcsec square subregions within the source regions shown in Figure 12, green 
and red referring respectively to the warm and hot gas components. For HH7, the squares 
apply to the HH7 bow shock, while crosses apply to warm gas ~ 20 arcsec behind the bow 
shock (HH8 and nearby warm gas). The black curve shows the expected behavior for gas 
with an initial ortho-to-para ratio of 0.25 (lower panel) or 0.4 (upper panel) that has been 
warm for a period of 150/n2(H) yr, where 100n2(H) is the density of atomic hydrogen. The 
orange curve shows the ortho-to-para ratio in LTE. The cyan curves and asterisks show the 
model predictions of Wilgenbus et al. (2000; see text). 



-43- 



5.8 micron 9.2 micron 




15 10 5 -5 -10 -15 15 10 5 -5 -10 -15 



Fig. 14. — Continuum maps obtained for several continuum wavelengths toward HH7-11 
Field 1, which contains the strong continuum source SVS13 at its center (i.e. at position 
(0,0)). Logarithmic contours are spaced by a factor 2^/^, with the second-highest contour 
corresponding to an intensity one-half the maximum. To the extent that SVS13 can be 
regarded as an unresolved point source at these wavelengths, these maps represent the beam 
profile resulting at the end of the entire data reduction procedure. 
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Fig. 15. — Comparison of the SVS13 spectrum obtained vi^ith Spitzer (heavy black hne) with 
that obtained with ISO (red, blue, green and magenta lines, different colors being used to 
show the different SWS orders). The S'j^zfeer spectrum obtained with the Short-Low module 
appears in the left panel, and that obtained with the Short-High and Long-High modules 
appears in the right panel. 
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Fig. 16. — Comparison of the H2 S(5)/S(3) line ratio and the total H2, H2O and CO surface 
luminosities (in erg cm~^ s~^) derived for the simple slab model (red; see text) with those 
obtained in the detailed calculations of KN96 (blue). 



